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Abstract. We report on absolute magneto-transmission experiments on highly- 
doped quasi-free-standing epitaxial graphene targeting the classical-to-quantum 
crossover of the cyclotron resonance. This study allows us to directly extract the 
carrier density and also other relevant quantities such as the quasiparticle velocity 
and the Drude weight, which is precisely measured from the strength of the cyclotron 
resonance. We find that the Drude weight is renormalized with respect to its non- 
interacting (or random-phase-approximation) value and that the renormalization is 
tied to the quasiparticle velocity enhancement. This finding is in agreement with 
recent theoretical predictions, which attribute the renormalization of the Drude weight 
in graphene to the interplay between broken Galilean invariance and electron-electron 
interactions. 
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1. Introduction 

The low- frequency dynamical conductivity, cr{(ju)^ a fundamental property of electrical 
conductors, is also central for the physics of graphene and for possible applications 
of this material in opto-electronics [H El El Hj. A tempting description of a{uj) in 
graphene refers to a simple model of noninteracting Dirac-fermion quasiparticles with 
the characteristic linear dispersion E(k) = ±^F^|k|, which at zero temperature occupy 
the electronic states up to the Fermi energy |£'f| = h.VY^/7T\n\^ where vy stands for the 
(Fermi) velocity parameter and n is the free-carrier concentration. Such an effective 
single-particle (ESP) picture correctly reproduces the measured strength of interband- 
absorption processes with the characteristic onset at 2|£'f| and the dispersionless 
universal amplitude Winter = (^uni = ^ I (4/i) [3 [6l [71 El [9l [IHl [H] • Further simple reasoning 
in terms of Drude approach yields the well-known expression for the integrated strength 
V = crintTa{^)d(jJ of intraband processes (so-called Drude weight): 



in which rUc = \Ey\/vy accounts for the effective cyclotron mass of carriers at the Fermi 
energy. 

The experimental measurements available so far [T2l [T3] are in conflict with this 
simple prediction. They suggest a suppression of the Drude weight up to 40 % (notably, 
in graphenes with high carrier concentrations), and therefore raise the issue of the 
applicability of the ESP picture to describe the dynamical conductivity of graphene. 
These doubts are additionally fostered by a startling report on non- vanishing absorption 
below the 2|£'f| threshold for interband processes [llj. 

Formally, the Drude approach follows from the commonly used random phase 
approximation (RPA) [14], which indeed grants Eq. ([T]), however, with the velocity 
parameter fixed at its bare value [15j v^^o = 0.85 x 10^ m/s - the value of 'Uf,o can 
be accurately estimated, for example, from calculations based on density-functional 
theory at the LDA level [161 [El IlH]. Markedly, v^^o is smaller than the apparent 
(renormalized) Fermi velocity, = (1.0 — 1.1) x 10^ m/s, which is derived from most of 
spectroscopy studies on graphene-based systems [HEOliaiESlESllHEaEni^ (and 
from beyond-LDA theoretical approaches). In fact, the Drude weight determined by the 
bare velocity, Vq = {vf,o/'^f)T^^ could partially explain the experimental claims [T2l [T3] 
of a suppressed V. However, the RPA approach leading to Eq. ([T]) with ^f replaced by 
Vy^o misses some important physics [15]. A more adequate formalism, based on many- 
body diagrammatic perturbation theory [15j, provides a more complex scenario in which 
self-energy (i.e., velocity enhancement) and vertex (i.e., excitonic) corrections compete 
with each other. Nevertheless, when the carrier concentration is sufficiently high, the 
approximate validity of the ESP model of noninteracting Dirac-fermion quasiparticles 
is expected to be recovered and Eq. ([T]), with the renormalized (apparent) velocity 
parameter, may be justified. Thus the issue of the strength of the Drude weight in 
graphene, one of its very basic electronic properties, remains a matter of controversy. 
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In this paper we report on testing the ESP model against the results of 
cyclotron-resonance (CR) absorption experiments, performed on highly-doped graphene 
specimens. Consistently following the ESP picture, we notice that the Drude weight 
is fully transformed into the strength of the CR absorption in the classical limit 
and, with a good precision, into the strength of inter-Landau-level transitions in the 
quantum limit, i.e., into a quantity which is accurately determined from our absolute 
magneto-transmission measurements. The carrier concentration is deduced from the 
analysis of the oscillatory broadening of the CR response with Landau-level filling factor, 
whereas the effective mass (and, in turn, the apparent Fermi velocity) is directly read 
from the dispersion of the CR with magnetic field. The Drude weight, as measured 
from the strength of the CR absorption, is found to be consistent with the estimate 
given by Eq. ([T]) in which the velocity parameter has the meaning of its apparent 
(spectroscopically measured or theoretically renormalized) value. 

2. Experimental results 

The investigated graphene sheet was prepared by thermal decomposition of the Si- 
terminated substrate of 6H-SiC. Subsequent hydrogenation of the "zero layer" yields 
a quasi- freestanding epitaxial graphene monolayer [29]. Such specimens are typically 
highly p-doped with the Fermi energy ^ —0.3 eV and carrier mobility /x ^ 
(2 — 3) X 10^ cm^/ (V.s) [30]. In this paper, we present data from one particular specimen. 




Photon energy (meV) 

Figure 1. The real part of the longitudinal conductivity axx{^^ B) (in units of (Juni) as 
extracted from the absolute magneto-transmission T{B) using the procedure described 
in Ref. [28]. Individual curves are plotted every 1 T. Thick curves are every 4 T. 
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Figure 2. Upper panel: Landau level fan chart with schematically-shown cyclotron 
resonance transitions in the quantum regime. Lower panel: A color plot of the real 
part of the experimentally determined longitudinal optical conductivity (Jxxi^^B). 
The dashed lines correspond to transitions between adjacent (hole) Landau levels in 
graphene, L-rn ^-m+i- In both panels, a Fermi velocity vp = 0.99 x 10^ m/s was 
considered. 



however, consistent results have been obtained on another similar sample. 

To measure the far-infrared transmittance, a macroscopic area of the sample 
4 mm^) was exposed to the radiation of a globar, which was analyzed by a Fourier 
transform spectrometer and delivered to the sample via light-pipe optics. The light 
was detected by a composite bolometer placed directly below the sample, kept at 
T = 1.8 K. From the substrate-normalized transmission T[B)^ the real part of the 
optical conductivity, [a^xi^^B)]^ was obtained for every separate magnetic field. 
The extraction of the conductivity from the transmission is described in detail in the 
Appendix. 
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Experimental data are presented in Figs. [T] and [2] in the form of the conductivity 
spectra and as a false-color plot of 3?e [axx{(^^ B)]^ respectively. In agreement with 
previous studies on highly-doped graphene [3ll [30], the magneto-optical response 
corresponds to the quasi- classical regime of the CR [32l [33l [34|, since the condition 
fiB ^ 1 applies. The CR line is approximately a single Lorentzian, its position in the 
spectra increasing linearly with 5, see Fig. [2l At low magnetic fields (in particular 
below 5 = 10 T), a departure from linearity is visible, see Fig. [3](a). In this regime the 
optical response is dominated by plasmons confined on the scale of a few microns due to 
graphene grain boundaries and/or substrate terraces. This regime has been discussed 
in detail elsewhere [35]. 

In high magnetic fields, the CR absorption still follows an overall linear-in-5 
evolution, consistently with moderate mobility of the sample. However, at the same 
time, it shows first clear indications of the approaching quantum regime, in which CR 
absorption corresponds to transitions between well-defined adjacent Landau levels (LLs). 
A number of CR studies of graphene in this regime has already been performed, see, 
e.g., Refs. [l9llMllZIll2ai23l[24l[25l[^ A typical response in the quantum regime 
closely refiects the properties of the LL spectrum: 

Ei = 8ign{i)vF^/2eh\i\B, i = 0, ±1, ±2 . . . . (2) 

Inter-LL excitations follow a clear \/S-dependence, and the CR absorption gains a multi- 
mode character [36] due to the non-equidistant level spacing. In the low-temperature 
limit we always get at least two CR absorption modes unless the Fermi level is placed 
precisely between two LLs. 

A closer inspection of the data in Figs. [Hand [2](a) reveals broadening of the CR line 
that appears at magnetic fields of 5^ ^ 27, 20, 16, 13.5 T and possibly also at 11.5 T. 
This broadening is best manifested as a decrease in the peak intensity, see Figs. [H or 
[21 whereas the total peak area remains fairly constant, see Fig. [3](b). This sequence 
of field values matches the simple rule mBm = const (with m = 3, 4, 5, 6 and 7) and 
allows us to identify the filling factors u = —4m, which correspond to the half-filled 
last occupied LL. The broadening at fields Bj^ thus refiects two CR modes, namely the 
transitions ^-m+i and Iz-m, with nearly equal intensities. It is worth 

noting that the CR width oscillating with the filling factor might also be reminiscent 
of another phenomenon - the diflFerence in the screening properties of a full or half- 
filled last-occupied LL in a 2D parabolic-band electron gas [37l [38l [39l [40]. Apart from 
the discussed weak oscillations in z/, the CR width does not exhibit any clear field- 
dependence, which would facilitate the identification of the main source of scattering in 
the sample, based on predictions of available theories, see, e.g., Refs. [HJ [42l [43] . 

The sequence of the CR line broadening = const is compatible with two 

possible scenarios. Either the Fermi level in graphene is pinned by electronic states 
in the SiC substrate, i.e., Ey = const = —v^V^eWmB^, or the carrier density n 
simply remains constant with varying B. In the latter case, the density is given by 
|n| = const = Nfm({B), where C(^) is the LL degeneracy eB/h and A^f = 4 is the 
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Figure 3. Part (a): Maximum of the CR absorption as a function of the magnetic 
field. The dashed hue denotes the theoretical position of CR for the effective mass 
of m = 0.058 mo- The departure of the CR position from the linear-in-5 behavior 
at low fields is consistent with the appearance of confined magneto-plasmons [35l |44] . 
Part (b): Drude weight (area) extracted using a Drude-Lorentzian formula (see Refs. 
[35l|44]) and the directly read-out amplitude of 5Re [axx{^^ B)]. The horizontal dotted 
lines show the expected ESP and "bare" Drude weights for a specimen with a 95% 
graphene coverage (see the main text), taking account the apparent (renormalized) 
vy = 0.99 X 10^ m/s and bare vy,o = 0.85 x 10^ m/s Fermi velocities, respectively. 



number of fermion flavors in graphene (i.e., a degeneracy factor associated with spin 
and valley degrees of freedom) . Even though charge transfer between the substrate and 
graphene cannot be excluded {Ey = const), this effect is negligible as testified by the CR 
weight being nearly constant with 5, see Fig. [3](b). Anyway, in both cases, the carrier 
density at 5 = and also at each 5^ is |n| = AemBj^/h = (7.9 ± 0.2) x 10-^^ cm~^. 
Importantly, this value of carrier density is extracted entirely from the LL occupation, 
i.e., independently of and E'p. This represents a significant advantage over previous 
studies of graphene in the quasi-classical regime [3ll [30]. We also note that the 
observation of individual filling factors in data, revealed as well-defined steps in position 
of the CR line in Fig. [21 serves as an indication of high homogeneity of the carrier density 
on the macroscopic scale. 

Using the ESP expression for the cyclotron frequency, Uc = eB/rUc^ we can extract 
the cyclotron mass at the Fermi level from the slope of the CR line directly: we find 
rric = (0.058 ± 0.01)me, where ttIq is the electron mass in vacuum. Using the "Einstein 
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relation", |£'f| = rricV^^ together with |£'f| = hv^ ^/7^\n\^ we obtain and E'p. We 
find Vy = (0.99 ± 0.02) x 10^ m/s, which is sensibly larger than the "bare" value 
^F,0 5 and |£'f| = (325 ± 5) meV which is in excellent agreement with the position 
of onset of the interband absorption (not discussed in this paper), which provides 
|£'f| = (320 ± 10) meV. According to the ESP model ([I]), the deduced Fermi level 
implies the expected Drude weight Vh/a^^i = 2\Ey\ = 650 meV. 

On the other hand, the Drude weight can be directly extracted from our data 
by estimating the area under [axx{(ju^ B)]^ which bypasses the use of Eq. ([T]). To 
account for the partial transfer of the Drude weight into the lower magneto-plasmon 
branch [45l l46l l47| . which is significant at lower magnetic fields, the Drude weight has 
been extracted by fitting with a standard Drude-Lorentz formula, see Refs. [35l [44] . 
assuming the confined-plasmon frequency of cjq = 6.5 meV [35]. This approach provides 
at higher fields results identical to the fitting with a simple Lorentzian curve. From this 
fitting, we obtained the Drude weight close to Vh/a^^i ^ 600 meV, see Fig. [3](b). 

To compare both estimates of the Drude weight, the following corrections should 
be taken into account: i) the coverage of the substrate by graphene is not full (see, 
e.g.^ AFM measurements on similarly-prepared samples [35j), ii) the hydrogenization 
process is not always complete - a small part of the sample remains covered only by 
the non-graphene "zero layer" [48j, iii) bilayer graphene may also appear at selected 
locations [49] and iv) the area below [axx{(^^ B)]^ extracted using the simplified 
Eq. IA.41 is suppressed with respect to the Drude weight by a few percent, see Appendix. 
All these corrections tend to decrease the Drude weight deduced from the area under 
[<^xx(^5 B)] by approximately 5-10%. To sum up, we find that our two independent 
estimates of the Drude weight agree with each other with a precision better than 10% 
and conclude that we do not observe any significant deviation from the validity of the 
ESP model ([H). Our findings thus do not support recent transmission studies [12], [13] 
in which a significant suppression of the Drude weight (in comparison with the ESP 
expectation) has been reported. The reason for this discrepancy remains unclear at the 
moment, although we believe that (at least a part of) the suppression of the Drude 
weight found in Ref. |12j might stem from the used normalization procedure, in which 
the Drude-type absorption of graphene in the regime of electron-hole puddles (when the 
sample is neutral on average) is neglected. 

3. Discussion 

The validity of the ESP model ([H) for the description of the strength of intraband 
absorption reminds us of Landau's theory of normal Fermi liquids, which establishes a 
"mapping" between a system of interacting fermions and a gas of weakly-interacting 
quasiparticles with renormalized parameters [2]. The renormalized quasiparticle 
velocity is precisely one of these parameters, which is completely controlled by the 
real part of the quasiparticle self-energy evaluated at the Fermi level E'p. Quasiparticles, 
however, do interact among each other and their interactions, which are described by 
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the dimensionless Landau parameters Fj'^, modify physical observables such as the 
compressibihty and the spin susceptibihty, or, more in general, the macroscopic response 
of the interacting electron liquid to external fields [14]. Our experimental test of the 
validity of the ESP model tells us that interactions among quasiparticles in graphene, 
which diagrammatically are encoded in "vertex" or "excitonic" corrections, are weak, 
at least at large carrier densities. 

Our experimental findings and the approximate validity of Eq. ([T]) can also be 
interpreted in terms of an enhancement of the Drude weight V with respect to its non- 
interacting value Vq^ which is found to be given by V/Vq ^ Vy/vy^o 1.2. This result 
is in excellent agreement with the theoretical predictions by Abedinpour et al. [15j. 
By taking into account screening at the Thomas-Fermi level, the authors of Ref. [15j 
predicted T> /T)^ ^ 1.25 at a carrier density |n| = 8 x 10''^^ cm~^ and for a dimensionless 
coupling constant o^ee ^ 0.5 - for definition of o^ee see, e.g., Ref. [Ij. Note that this is 
the appropriate value of o^ee for our sample, as deduced from a recent ARPES study 
carried out by Bostwick et al |50j in a nearly identical sample. We stress that the 
theoretical proof of the approximate validity of Eq. ([T]) given in Ref. [15j is far from 
trivial: it is based on a fully-microscopic diagrammatic approach in which substantial 
cancelations between the self-energy and vertex corrections play an essential role. In 
particular, the authors of Ref. [I5] demonstrated the smallness of excitonic corrections 
to the intraband response of a doped graphene sheet, which tend to suppress the Drude 
weight - see Eqs. (21) and (23) in Ref. p!5]. Our experimental study corroborates this 
finding too. 

For the sake of completeness, we should note that the Drude weight V is strictly 
fully transferred into the CR strength only in the quasi-classical limit. In the quantum 
regime, the total CR weight Vcr may differ from V. When the Fermi level is placed 
just in between two LLs = 4m + 2, m = 1, 2, 3 . . .), we still have V = Pgr, but for 
the half-filled last occupied LL (|z/| = 4m, m = 1, 2, 3 . . .) we obtain Vqr < T> fST] : 



However, one can easily verify that this suppression, which is related to the transfer of 
the oscillator strength between intra- and inter-band absorptions, remains very small 
for m > 2 (below 2%). The approximation V = Vqr is thus very well applicable as long 
as the zero LL is not involved in the CR absorption. 

4. Conclusions 

To conclude, an absolute magneto-transmission experiment on highly-doped quasi- 
freestanding graphene allowed us to follow the oscillator strength of the cyclotron 
resonance, which exactly matches the zero-field Drude weight. We found that this 
weight remains very close to V = 2|£'F|cruni/^ = e'^vv a/^I^I/(2^) = 7r|n|e^/(2mc) and 
that a single-particle-like picture is fully sufficient to account for our data, provided 
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that the renormahzed (quasiparticle) parameters E'p, '^f and/or rric are considered. This 
imphes that the Drude weight in a doped graphene sheet is enhanced with respect to 
its bare value due to electron-electron interactions and that the enhancement is tied to 
the Fermi velocity enhancement [52l [53]. Our findings are in excellent agreement with 
recent theoretical work p!5], which predicted an enhancement of the Drude weight in 
doped graphene stemming from the interplay between broken Galilean invariance and 
electron-electron interactions. 
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Appendix A. Extraction of the optical conductivity from transmission 
spectra 

The experimentally obtained substrate-normalized transmission T{B) spectra are used 
to extract the optical conductivity Re{axx{(jU^ S)}, of the graphene on the SiC substrate 
at every field. Figure lAll shows the complex transmission and (internal) refiection 
coefficients for our experimental "vacuum (v) / graphene (g) / substrate (s) / vacuum 
('u)" geometry which are: 

2 . 2Ns . 2 



f = f = — f = 



^sv = T7~r~7' '^sgv = J , ^ , ry ^ = exp{iu Ngdg / c) . (A.l) 



N,-l _ N,-l-Zoa 

Here ^5 = 77,5 + ikg is the experimentally known complex refractive index and dg is 
the thickness of the substrate, Zq = 377 Q is the impedance of vacuum and a is the 
two-dimensional optical conductivity of graphene. 

In the derivation of the ratio of the transmission coefficient of the sample to the one 
of the bare substrate all internal reffections in the substrate must be taken into account. 
The experimental spectral resolution was reduced to 4 cm~^ in order to suppress the 
Fabry-Perot interference in the SiC substrate. In this case, the internal reffected rays 
add incoherently. By plane-wave counting, the experimental transmission coefficients of 
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Figure Al. The schematic representation of the sample and the definitions of the 
complex coefficients of Eqn. (jA.ip . A near normal incidence was used, the rays on the 
figure are inclined for clarity. 



the bare substrate and the sample are derived in terms of the complex coefficients: 

|2 u ^ |2 



T = 

vsv 



- T 

2|2' -^vgsv 



\tvgstsv^\ 



|2 • 



1 — r r 1 
The substrate-normahzed transmission foUows from Eqns. (lA.ip and (IA.2p : 



(A.2) 



T 



T 

-J- 1 



vgsv 



T 

-L- n. 



\Ns + l 


4-|iv,-ir 




4 


\{N, + l)(iV, + 1 + Zo(t)|2 -\{N,- 1){N, - 1 - Zo(t)|W 



In the experimental spectral range, kg <^ rig^ therefore kg is neglected in all complex 
coefficients f lA.ll ). except in r. We introduce the absorption coefficient of the substrate 
a = |rp and the substrate- vacuum amplitude reflectivity r = {ug — l)/(^s + 1) in Eqn. 
( 1A.3I ) and after some algebra the relation between the substrate-normalized transmission 
and the complex optical conductivity is found: 

-1 



T = 



1 + 



1 



a'^r^ 



Re{Zoa} + 



1 



1 



2 2 



Q,2^4 



Zo(t\ 



. (A.4) 



+ ""^ ' (n, + 1)2 1 

Evidently, from only the experimental transmission spectra the complex conductivity a 
can not be extracted. However, note that the prefactor in the term depending on Re{cr} 
is significantly larger than the prefactor in the term depending on |crp, therefore the 
transmission is determined predominantly by Re{a}. In this case replacing |ap with 
Re{cr}^ is a reasonable approximation. In magnetic field, Eq. f lA.4p remains valid for 
each circular polarization separately, if a is substituted with a±. For unpolarized light 
used in our experiment, the transmission is an average of the two. In order to extract 
the diagonal conductivity at finite field, we made a further simplification by ignoring 
the relatively small effect of on the transmission and substituting a with Re{axx} in 
Eqn. f lA.4p . To estimate the error introduced by these approximations, we have tested 
the simplified version of equation IA.4I using the ideal classical CR absorption described 
by the standard analytical formulae for and axy. We have found that the area under 
Re{axx}^ extracted using the above described procedure, is suppressed with respect to 
the Drude weight by roughly 5%. 
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